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Calculation of two-loop self-energies in the electroweak Standard Model* 
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Motivated by the results of the electroweak precision experiments, studies of two-loop self-energy Feynman 
diagrams are performed. An algebraic method for the reduction of all two-loop self-energies to a set of standard 
scalar integrals is presented. The gauge dependence of the self-energies is discussed and an extension of the 
pinch technique to the two- loop level is worked out. It is shown to yield a special case of the background-field 
method which provides a general framework for deriving Green functions with desirable theoretical properties. 
The massive scalar integrals of self-energy type are expressed in terms of generalized multivariable hypergeometric 
functions. The imaginary parts of these integrals yield complete elliptic integrals. Finally, one-dimensional integral 
representations with elementary integrands are derived which are well suited for numerical evaluation. 



1. INTRODUCTION 

The beautiful results of the LEP experi- 
ments have shown that the electroweak theory 
has a predictive power like QED several decades 
ago. It is to be expected that eventually the elec- 
troweak theory will provide high precision pre- 
dictions for many experiments in the present and 
near future. One may in particular think in this 
respect of the measurement of the Z mass and 
width, a better determination of the W mass and 
a direct indication of the top mass. 

These developments will require in the future 
two-loop calculations in the electroweak theory. 
In the study of the two-loop contributions the 
self-energy diagrams play a central role. The 
present paper reviews the problems one encoun- 
ters when considering self-energies in the Stan- 
dard Model (SM) and discusses recently devel- 
oped methods to tackle them. 

The choice to treat first the self-energies is a 
logical one. They are the simplest two-loop di- 
agrams which yield a universal contribution to 
all two-loop processes. Moreover, from one-loop 
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calculations we know their importance. Several 
results for two-loop self-energies are known, in- 
volving however approximations. They concern 
the limiting cases of a heavy fermion doublet [ 2 1 , 
a heavy top quark |^ and a large Higgs mass [ i | . 
However, in the electroweak theory also the W 
and Z bosons have a non-negligible mass, and 
one is in general faced with two-loop self-energies 
containing massive propagators. 

Let us summarize the problems which arise 
when evaluating two-loop self-energies. The first 
problem is the plethora of Feynman diagrams 
contributing, the second one the evaluation of 
tensor integrals in terms of scalar integrals and 
the third one the calculation of the scalar inte- 
grals themselves. When discussing the correc- 
tions due to two-loop self-energies, one further- 
more has to deal with the problem that these 
contributions are in general gauge dependent. 

In order to handle the large number of Feyn- 
man diagrams needed for the evaluation of two- 
loop self-energies, an algebraic approach is chosen 
allowing for a high degree of automatization. It 
involves a method for the decomposition of tensor 
integrals. Whereas for one- loop diagrams this is 
a well-known procedure it was only recently 
developed for two-loop self-energies ^. In this 
way all two-loop self-energies can be reduced to 
a set of standard scalar integrals. The task of 
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evaluating any two-loop self-energy involving in 
general several thousands of Feynman amplitudes 
is therefore reduced to the problem of calculat- 
ing four different types of two-loop scalar inte- 
grals. As an example, we treat light fermion con- 
tributions to the self-energy of the .Z-boson. The 
results given in a minimal basis of standard in- 
tegrals allow to study the gauge dependence of 
the considered quantities directly at the algebraic 
level. 

At one-loop order the self-energies are fre- 
quently used as building blocks to define running 
couplings or to parametrize electroweak radiative 
corrections. The need for self-energies as building 
blocks with suitable theoretical properties seems 
to be even more important at the two-loop level, 
where so far no calculation of a complete pro- 
cess exists. In the one-loop applications it was 
found that the self-energies evaluated in the class 
of i?5-gauges are not adequate as building blocks 
due to their gauge dependence and unsatisfactory 
high energy and UV behavior. Many proposals to 
modify these self-energies aimed on eliminating 
their gauge parameter dependence. In particu- 
lar, the pinch technique (PT) was found to 
yield results with decent properties. Recently it 
was shown ||^ that the background-field method 
(BFM) j|] offers a wider framework in which 
Green functions possessing desirable theoretical 
properties are directly derived from a gauge in- 
variant effective action. In this paper we work 
out an extension of the PT to the two-loop level 
and show that it corresponds to a special case of 
the BFM results. 

The evaluation of the two-loop scalar integrals 
is in general needed for non-negligible masses 
and an arbitrary dimension D. When expanding 
around S — {4 — D)/2 ^ 0, one gets in these cases 
not anymore results in terms of (poly) logarithms. 
To perform nevertheless an evaluation of the mas- 
sive two-loop diagrams, in essence two approaches 
are followed: expansions for small and large exter- 
nal momenta and numerical integration. In both 
strategies this paper presents recent new results. 
On one hand, the use of multi-variable general- 
ized hypergeometric functions, i.e. multiple se- 
ries, is a new development. On the other hand, 
one-dimensional integral representations are de- 



rived, which are a good alternative to the exist- 
ing two-dimensional integrals. It turns out that 
the derivation of one-dimensional integrals is pos- 
sible for all two-loop diagrams which contain a 
one-loop self-energy insertion. 

For the imaginary parts of the scalar integrals 
simple analytic results in terms of complete ellip- 
tic integrals are derived. 

The outline of the paper is as follows. 

The algebraic method for reducing two-loop 
self-energies to standard scalar integrals is de- 
scribed in section 2, whereas in section 3 the ques- 
tion of gauge invariance in the framework of PT 
and BFM is discussed. The next section deals 
with an analytic approach to scalar self-energy 
integrals leading to generalized hypergeometric 
functions. Section 5 focuses on the imaginary 
part of the scalar integrals, which are then ex- 
pressed in terms of elliptic integrals. The one- 
dimensional integral representations are derived 
in section 6. In section 7 we give our conclusions 
and an outlook. 

2. ALGEBRAIC CALCULATION OF 
TWO-LOOP SELF-ENERGIES 

2.1. Classification 

As the first step in the evaluation of two- 
loop self-energies, the relevant Feynman diagrams 
have to be generated. These can be classified ac- 
cording to their topologies. In fig. |l| some exam- 
ples of irreducible two-loop self-energy topologies 
are shown. A complete set including also tadpoles 
and topologies which factorize into one-loop con- 
tributions is given in . Inserting fields into the 
topologies in all possible ways and applying the 
Feynman rules yields all relevant Feynman ampli- 
tudes. The generation of the Feynman diagrams 
and amplitudes is carried out automatically with 
the program FeynArts pO| . 

In the following we will use the shorthand no- 
tations 

(( ))=/ ^ / ^ ( ) 

where qi and q2 are the integration momenta of 
the loop integrals and /i is an arbitrary reference 
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Figure 1. One-particle irreducible topologies of two-loop self-energies. 



mass, and 
dij...i = 



Here ki is the momentum of the i-th propaga- 
tor and rrii its mass. In the most general case 
corresponding to the first topology in fig. |^ five 
different propagators appear. Their momenta are 
related to the integration momenta qi and q2 and 
the external momentum p via 

ki = 9l, ^2 = 91 +P, ^3 = 92 - 91, ^4 = 92, 
^5 = 92 +P- (1) 

A general Feynman amplitude takes the form 

{{f{h,m,,D)di2...i)), (2) 

where f{ki,mi,D) is an expression depending 
on the momenta, the particle masses and the 
space-time dimension D and containing the whole 
Lorentz and Dirac structure. 

As will be discussed in the next section, all 
two-loop self-energies can be reduced to a class of 
standard scalar integrals which we call T-integrals 



rp / 2 2 2 

= {{■ 



kl - m 



», (3) 



where the denominator is of the same form as in 
the Feynman amplitude. The double index no- 
tation Tiji2 is used here to indicate that the 
indices of the T-integrals refer to the correspond- 
ing momenta fci ^ , fci2 , . . . ^i, . In the following the 
momentum and the masses will only be explic- 
itly written as arguments if confusion is possible. 

The topologies depicted in fig. ^ correspond to 
the scalar integrals T'12345, ^11234, ^1234, 1234 and 
7^1134, respectively. The analytical expression for 
T'11234 can be obtained from T1234 by partial frac- 
tioning or taking the derivative with respect to 



nil ■ Other integrals with higher powers of propa- 
gators are treated in the same way. For the gen- 
eral case one therefore needs to evaluate only four 
different types of two-loop scalar integrals. 

2.2. Reduction to standard scalar integrals 

It is convenient to begin with a decomposition 
into Lorentz scalars. For the gauge boson self- 
energies it reads 



Pp.Pv 



from which the transverse part J^t{p'^) and the 
longitudinal part S]i(p^) can be obtained 



1 



D - 1 



^tiu{p)- 



For all other types of self-energies scalar quanti- 
ties can be extracted in a similar way ^J. 

The contraction of Lorentz indices, reduction of 
the Dirac algebra and evaluation of Dirac traces 
can be worked out like in the one- loop case, 
e.g. with the program FeynCalc [Tl| ]. Since we 
deal with scalar quantities, this results in scalar 
products of momenta (fci'fcj), [ki-p), in the nu- 
merator of the Feynman amplitude. The denom- 
inator is unchanged. We now implicitly use mo- 
mentum conservation and express all scalar prod- 
ucts as sums of momentum squares, e.g. 



{kl 



kl-p^ 



(4) 



Subsequently all kf appearing both in the numer- 
ator and the denominator are canceled via 



(5) 



In general there remain kf in the numerator 
which cannot be canceled, e.g. for 



I={{kldi2ZA)). 



(6) 
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We proceed in this case by writing 

kl = {kl - ml) + (ml+p^) + 2{p ■ k^). 



(7) 



The first term can be canceled with the appropri- 
ate propagator factor in the denominator. The 
second term contains no integration momenta. 
We therefore focus on the integral 



(8) 



and perform a tensor decomposition for 5']'234- 

In contrast to the one-loop case the decompo- 
sition with respect to the external momentum, 
i.e. the ansatz Istl 



(9) 



does not lead to simpler integrals here. In order 
to determine the scalar quantity S{p'^), one has 
to contract with yielding terms in the numer- 
ator which cannot be canceled. This fact is due 
to the structure of the topology associated with 
•^1234 third topology in fig. ^ which contains 
a four- vertex with three inner lines. This is a typ- 
ical feature of two-loop topologies not present at 
the one-loop level. 

Instead, we perform a decomposition with re- 
spect to a subloop. We write 



{k^d34/ 



— A:^s(fc]^), 



(92 - ^i)^ - "^i 



91 



(10) 



where the last equality follows from the fact that 
the tensor integral depends only on fc^. The 
quantity s(fcf ) is given by 



s{kl) 



1 



2 {{kl ■ ki)d3i) 



{{kl ■ kii)d 



1'34/ 



(11) 



The factor l/kf can be written as a massless prop- 
agator, which we have indicated with a prime at 
the corresponding index in (ii'34. Inserting this 
into (||) yields 

Pm'5'i234 = {{{P ■ kl){ki ■ fc4)dl'1234))- (12) 

The scalar products in the numerator of this in- 
tegral can now be expressed in terms of squared 



momenta which can be canceled. Therefore all 
momentum dependent terms have been removed 
from the numerator of the Feynman integral while 
its denominator has retained its original struc- 
ture. It can further be simplified by performing a 
partial fractioning for the propagator factors dyi 
carrying the same momentum. 

It can be shown |^ that in a similar way all pos- 
sible integrals appearing in calculations of gen- 
eral two-loop self-energies can be reduced to T- 
integrals. 

The result of the algebraic calculation in gen- 
eral contains several scalar integrals depending on 
different arguments. These are not necessarily in- 
dependent of each other, i.e. there exist relations 



c"r" = 



(13) 



where , . . . , c" are polynomials in p^ . We use 
these relations together with the symmetry prop- 
erties of the T-integrals to eliminate as many in- 
tegrals as possible. This leads to a minimal set of 
scalar integrals which are algebraically indepen- 
dent of each other. The result expressed in this 
minimal basis is very transparent and directly dis- 
plays certain properties of the considered quan- 
tity. 

For example, Slavnov- Taylor identities can be 
checked at a purely algebraic level without having 
to use any explicit expression for the T-integrals. 
When adding up the results for the relevant am- 
plitudes, the coefficient of every standard integral 
exactly adds up to zero. This has been demon- 
strated in for Slavnov- Taylor identities involv- 
ing several thousands of Feynman amplitudes. 

The algorithms outlined here have been imple- 
mented into the program TwoCalc [l^ ] which car- 
ries out the algebraic calculation fully automati- 
cally. 

2.3. Results for the Standard Model 

Another feature which can directly be read off 
from the algebraic result is its gauge dependence. 
We work in a general i?^ -gauge specified by one 
gauge parameter {i — 7, Z, W) for each vec- 
tor boson and write the lowest-order gauge boson 
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Figure 2. Light fermion contributions to the two-loop self-energy of the Z-boson. 



propagators as 



i(l - l/^i)kfj,ki, 



fc2 



fc2 



.(14) 



As a simple example, we consider two sets of light 
fermion graphs contributing to the Z self-energy. 
They are shown in fig. |2| A complete listing of 
the light fermion contributions to the gauge boson 
self-energies was given in . All fermions except 
the top-quark can be treated as light fermions, 
i.e. their masses are small compared to those of 
the Z- and VF-bosons and can therefore be ne- 
glected. For definiteness, we choose d- and u- 
quark. 

We find for the transverse part of the three 
"neutral current" diagrams depicted in fig. |3a 



,ZZ,(2) 



^T,z '{p") = e''Cz[fziAo,Bo) 

+ ii-6D + D^)Tiy4,im%) 
+ {8-4D + D^)T23'4,ip^;m%) 
-{A-8D + D^){ml + p')ri'2'34' (p'; m|) 
4 



-D{m\ 



4m|p2 - m%p^D/2 + /) 



2'34'5' 



(15) 



where Cz is a dimensionless constant and 
fz{Ao,BQ) represents a function containing only 
scalar one-loop integrals. The result for each sin- 
gle diagram of fig. ||a depends on the gauge pa- 
rameter However, as is seen in (|l|), in the 
sum of these diagrams the gauge parameter has 
canceled, i.e. this contribution is gauge indepen- 
dent within the class of i?{-gauges. 

Fig. Hb represents a set of "charged current" 
diagrams where a is exchanged. We have 



only drawn one diagram of each type. The result 
for the altogether 15 amplitudes is given by 



fiAo,Bo) 



"'■W"''Z 

+g{Ao,Bo;^w) + F{T) 



(16) 



Here C is again a dimensionless constant, 
/{AqjBq) and g{AQ, Bq; represent terms 
which only involve scalar one-loop integrals, and 
F{T) and G{T; ^w) contain the scalar two-loop 
integrals. In contrast to (15), this result is gauge 
dependent. For illustration, we only give here the 
explicit form of the gauge parameter dependent 
function G{T;^w)- 



T23'4'{p^;m'^/^w) 
\t,w ' 



-2 



(3 - 2D)p'^ 



4 



(-^ - l)' - 2(-^ + 3 - 2D)m^^p 

\l;W ' ' 



2 , 



w) 



(17) 



The gauge parameter appears both in the argu- 
ments of the scalar integrals and in their coeffi- 
cients. The other terms contributing to ( |l6| ) were 
explicitly given in Ref. 

The characteristic feature of the result ( p^ ) is 
the factor [j? — rn^) multiplying the gauge de- 
pendent part G(T; ^w) which contains the generic 
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Figure 3. Two-loop vertex and box graphs containing propagator-like pinch parts. 



two-loop integrals. This reflects the fact that the 
pole position of the two-loop propagator is gauge 



independent. Note that for p — 



there is 



still a remaining gauge dependence in the func- 
tion g{Ao, Bo;^w)- It cancels with the gauge de- 
pendence of the terms generated by inserting one- 
loop countcrtcrms into the one- loop diagrams. 



GAUGE INVARIANCE 
LOOP SELF-ENERGIES 



OF TWO- 



In order to investigate physical effects due 
to two-loop contributions of propagator type, it 
would be desirable to arrange these contributions 
in such a way that they form building blocks 
with suitable theoretical properties. At one-loop 
order, the pinch technique (PT) was devel- 
oped to eliminate the gauge parameter depen- 
dence by shifting contributions between differ- 
ent Green functions. It was found that the new 
"Green functions" obtained in this way fulfill sim- 
ple Ward identities and in comparison to their 
i?j-gauge counterparts possess desirable proper- 
ties such as improved IR and UV properties and 
a decent high-energy behavior. 

However, it was recently shown that the 
requirement of gauge parameter independence 
is not the crucial criterion for obtaining well- 
behaved Green functions. The background-field 
method (BFM) provides a more general frame- 
work in which the Green functions are directly 
derived from a gauge invariant effective action. 
Their desirable theoretical properties are a conse- 
quence of gauge invariance and hold for all values 
of the quantum gauge parameter ^g. In ||^ the 
BFM was applied to QCD and the SM. It was 
shown that the BFM includes the PT results as 
the special case (,q = 1- Viewed from the frame- 
work of the BFM, the PT results therefore are 
not gauge independent but correspond to a cer- 




Figure 4. The pinch part of the first diagram in 
fig. 3. 

tain choice of the quantum gauge parameter ^g. 

In order to discuss this issue at the two-loop 
level, we present here for an example the applica- 
tion of both approaches. In contrast to the BFM, 
the PT has so far been restricted to the one-loop 
level only. We therefore work out an extension 
to the two-loop case. We treat the "charged cur- 
rent" light fermion contributions to the two-loop 
Z self-energy. The result in the i?{ gauge was 
given in (Iq). 



3.1. Pinch technique at the two-loop level 

In the PT, the gauge parameter dependence of 
self-energies is canceled by combining them with 
propagator-like pieces extracted from vertex and 
box diagrams which contribute to a gauge inde- 
pendent S-matrix element. The relevant types of 
two-loop vertex (containing also the wave func- 
tion corrections) and box diagrams for a four- 
fermion process are drawn in fig. |^. 

Parts of these vertex and box graphs in which 
the propagators of the external fermions have 
been canceled become propagator-like. Such a 
"pinch part" is depicted in fig. ^. It corresponds 
to the first diagram in fig. ^, in which the propa- 
gator associated with the final state fermions has 
been removed. The precise definition of the pinch 
parts a priori is not unique. In analogy to the one- 
loop case, we identify the pinch parts as terms 
emerging when a longitudinal momentum fc^ of 
an elementary vertex or propagator is contracted 



7 



with a coupling matrix 7'^ of the external fermion 
current. This gives rise to an elementary Ward 
identity, i.e. it can be written in terms of inverse 
fermion propagators and mass terms 



^m2) + mi-m2,(18) 



where pf is the momentum of the external 
fermion. The first term on the right-hand side 
cancels the fermion propagator and therefore 
gives rise to the pinch part, whereas on-shell the 
other terms yield contributions proportional to 
the fermion masses. 

The /c^-terms originate from the longitudinal 
part of the gauge boson propagators in the R^- 
gauge (see (|l4|)) and from the three-gauge-boson 
vertices. Up to group theoretical factors these 
vertices can be decomposed into 



P 
F 



where 
iq + k) 

'2qf_cgua - ^qiygf^a - (2k + q)agny (19) 



The momenta (5 + fc)^ and in F^^^ have the 
same form as those arising from the longitudi- 
nal part of the relevant gauge boson propaga- 
tors. Therefore T^^,^ also gives rise to pinch 
parts. Note that in the first diagram of fig. ^ only 
the part of F^^^^ directly acting on the external 
fermion line yields a pinch part. 

Extracting the pinch parts of the diagrams rep- 
resented by fig. 1^ and combining them with the 
usual self-energy we find a result of the form 



E 



(p^)-f4^)p,„(p^) = 
4(/(Ao, Bo) + .g(Ao, So; iw) + F{T)) . (20) 



Here S^^'*'^'' (p^) denotes the usual self-energy in 
the i?5-gauge as given in (|l6|), E^^p-^Jp^) are the 
pinch parts of the one-particle irreducible two- 
loop triangle and box graphs, /(Ao, Bq) and F{T) 
are gauge independent functions of scalar one- 
loop and two-loop integrals, respectively, whereas 
5(^0, -Bo; Civ) still depends on the gauge param- 
eter. This means that only the gauge parame- 
ter dependence of the terms containing the scalar 
two-loop integrals has disappeared, while the con- 
tributions involving the one-loop integrals remain 



gauge dependent. This feature was in fact antici- 
pated from the structure of the result (^6|). Since 
the pinch parts extracted from the two-loop tri- 
angle and box graphs are by construction propor- 
tional to (p^ — TO^), they cannot cancel the gauge 



dependence still present in (^^ for p — ra^ ■ 

The complete two-loop contribution to the 
gauge boson propagators in the PT is given as 



(p2 



(21) 



(p^ — m^)^ 

where the first line contains the one-particle irre- 
ducible contributions of (|2^) and the second line 
the reducible ones. For the Z propagator one 
also has to take into account the mixing with the 
photon. Since E^pj„(p^) originates only from 
reducible graphs, it is not the full one- loop pinch 
part. In particular, it contains no contributions 
from one-loop box diagrams. However, these are 
necessary to cancel the gauge parameter depen- 
dence of the one-loop gauge boson self-energy, 
I]^«^(^)(p2) The reducible contribution in (1^) 
is therefore also gauge parameter dependent and 
does not correspond to the one-loop self-energy 
in the PT, S^^'^^\ 

To identify the quantity which has to be re- 
garded as the two-loop self-energy in the PT 
framework, we write the propagator (|l]) as 



PT,(2) 



This yields 



,PT,(2) 



ip' 



.(2) 
■'T.Pi 



.ip') 



TO' 



ip') 



^(1) 
-•T.Pin 



ip' 



(22) 



The two-loop self-energy constructed via the PT 
therefore involves both one-particle irreducible 
and reducible contributions. 

We have evaluated this quantity for the 
"charged current" light fermion contributions to 
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the Z self-energy and found that it is in fact in- 
dependent of the gauge parameters. The result 
will be given in the next section. 

3.2. The background-field method 

The BFM is a technique for quantizing 
gauge theories without losing explicit gauge in- 
variance. This is achieved by decomposing in 
the Lagrangian the usual gauge field into a back- 
ground field Y and a quantum field V . By adding 
a suitable non-linear gauge-fixing term, an effec- 
tive action r[y] is obtained which is invariant 
under gauge transformations of the background 
fields. While the background fields are associated 
with the tree propagators, the quantum fields ap- 
pear only inside loops. In r[y], only for the quan- 
tum fields a gauge fixing is needed. We specify it 
by the quantum gauge parameter ^q. 

The BFM is valid for all vertex functions to 
all orders of perturbation theory. Its approach is 
very different from the PT. While the PT seeks to 
obtain suitable "Green functions" by eliminating 
their dependence on the parameters of the i?^- 
gauge, the BFM is based on the gauge invariance 
inherent in r[T^]. The fact that the vertex func- 
tions fulfill simple Ward identities and possess 
other desirable features is a direct consequence 
of gauge invariance in the BFM |^ , whereas the 
PT provides no explanation for these properties. 

The comparison of these methods at the two- 
loop level clearly shows the power of the BFM 
approach. In the PT an involved rearrangement 
was necessary between different Green functions 
contributing to a specific two-loop process and 
between reducible and irreducible contributions. 
Contrarily, the vertex functions of the BFM are 
uniquely given in terms of Feynman rules. One 
simply has to calculate the "charged current" di- 
agrams shown in fig. ||b using the BFM Feynman 
rules. The result is just an ordinary self-energy 
which is guaranteed to be process independent. 
For the special case = I it reads 



4(1 - i?)/Bo(/; 0, 0)So(/; m^^,, m^,) 



H — 2^23' 4' {p ■,m^,) 
-2{3-2D)Ti23'4'{p'';ml,,', 



>^Tii23'A>{p ] 771^,7111^, m^)j ,(23) 

where Yj^'^^^\{p'^) is the conventional self-ener- 
gy in the 't Hooft-Feynman gauge, i.e. the result 
given in (^ with ^v/ = 1. The result (^ exactly 
coincides with the PT expression ( |2^ ) . Therefore 
we find here the same connection between the PT 
and the BFM as in one-loop order, i.e. the BFM 
contains the PT results as a special case. 

We would like to stress that the case = 1 
displayed here is not distinguished in the BFM. 
All properties derivable from the gauge invariant 
effective action, e.g. the simple BFM Ward iden- 
tities, hold for all finite values of ^g. As an ex- 
ample, we list a Ward identity valid for the Z 
self-energy in the BFM to all orders 



(24) 



'"W r 



c 13'4 



where the conventions of |^ are used. We checked 
the validity of this Ward identity at the two-loop 
level by explicit computation for arbitrary ^q. 
The Ward identity fulfilled by the conventional 
two- loop Z self-energy was given in It in- 
volves in addition the self-energy of the unphysi- 
cal scalar x a-nd reducible contributions. 

In analogy to the one-loop case, this means 
that the special case singled out in the PT is not 
unique. When studying applications of two-loop 
self-energies, it has to be investigated whether 
these are infiuenccd by the ambiguity correspond- 
ing to different choices of ^q. 

4. ANALYTIC APPROACHES AND HY- 
PERGEOMETRIC FUNCTIONS 

4.1. Introductory remarks 

As explained above, the algebraic reduction 
of two-loop self-energies leads to four types of 
two-loop scalar integrals, T234, 7x234, T12345 and 
T134. The last integral follows from T234 by set- 
ting p2 = (see below), so one needs to study 
three scalar two-loop integrals. For convenience, 
we will in the following relabel the integral r234 
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as Ti23(p^;m?,m2,TO|). 

In this section we discuss analytic results in 
an arbitrary number of dimensions D for these 
integrals. They will be in the form of generalized 
hypergeometric functions, that is multiple series 
of ratios of and the masses. We also discuss the 
expansion of these results around D — A. Only in 
those cases which fulfill the condition 



P 



(mi ±m2±TO3) ) =0(25) 



for the three-particle cuts, the finite part can be 
expressed in terms of polylogarithms. A typical 
example is given below: 

^ -25 

Ti23(p';0,0,m2) 



/ 9 

/ m 



r2(i-,5)r((5)r(-i+2<5) 



TO 



r(2-<5) 
1 1 

2^ ^ ^ 



2^^1(^,25-1; 2 -5; z) 



2 4 



3 



;C(2) 



13 , ^ 

- y2 + L«2(^) 



1 



- z log(l - z) 



(26) 



where z = /m?, L,n — j + ln(TO^/47r/x^), 
7 is the Euler constant and ({2) = 
Another example is the massive vacuum dia- 
gram Ti23(0; mf , TO2, m§) which in D dimensions 
is given in terms of four AppcU functions but in 
four dimensions contains only dilogarithms [p^ . 
In these specific cases one works with Feynnian 
parameters or dispersion relations [p^ . In the 
general mass case it turns out that x-space tech- 
niques and in particular the Mellin-Barnes repre- 
sentations are also useful. In the following subsec- 
tions we will give an example of the application 
of dispersion relations and of the Mellin-Barnes 
representation. 

4.2. The integral Tiaa 

Let us first start with the rederivation 
of the so-called London transport diagram 
Ti23(p^; TO^, m|, TO-g) by means of dispersion re- 
lations. The imaginary part or the discontinuity 
ATi23 is given by 



ImTi23{p'^;ml, ml, ml) 



ATi23{p^;ml,ml,ml) ^ 
2i 

-ni^np'r ~1 e (P' - (mi + m2 + m3)^) 



r2(2-2(5) 



ds 



Xi-\s,mlml)Xi-\p\s,ml) 



2U-<5 



= y- / dsABo{s:ml,ml)ABo{p'^:,s,mj) 
in J 

(m2+m3)2 

where A(a, b,c) = {a — b — c)^ — 46c is the Kallen 
function. The dispersion relation reads 

rri / 2 2 2 2\ 

Ti23(P ;m^,m^,m^) = 

oo 

If 1 

— / dz -ATi23{z;ml,mlml) 

ZTTl J Z — p'' 

(mi+m2+m3)^ 



J ds ABo{s;ml,ml) 



(m2+m3)2 

dz 



z — p^ 



ABo{z; s, ml). 



(28) 



(v^+mi)2 

Using the expansion 

1 I -r-^ f p'' 



r ( 7 



^ fc=0 ^ 



we perform first the integration over z 

„ oo 

dz 



(29) 



A = Y.(p')' f 



^k+2~S 



{z - u)'/^-\z - v)'/^-' 



= ^{p^)^u-^-^B{k + 5, 3/2 - 5) 

k=0 

X2Fi{S-l/2,k + S;k + 3/2;v/u) (30) 

where u = (mi-t-Y^)^ and u = (mi — -^5)^. One 
can transform the Gauss hypergeometric function 
2-^1 using relations which one can find in lead- 

4) 

v<5 



ing to ri23(p^; m^, m2, TOg"* 



(27) Ti23(p2;TO?,m2,TO2) = -(4V) 
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m-sns-i) 



r(2-2(5) 



ds 



(m2+m3)2 



2 / 

nil I 



-5 



F4(l,2-<5;2-5,2-5;^A) 
s s 



F4{1,6;2-S,S: 



p m 



(31) 



Since in ( |2q ) the z integration up to some factors 
represents the one-loop self-energy, we find as a 
byproduct 

" =T{S-1) 

-s 



9 / 



Fiil,S;2-S,5;- 



(32) 



Using the definition of the F4 functions we 
can easily perforin the integration over s. After 
some manipulations we obtain the result for the 
London transport diagram which agrees with the 
one derived using x-space technique and Mellin- 
Barnes representation jl^]. The result is pre- 



mf /m 



3 ' 



sented in table (|l|), where Zi 
1, 2 , Z3 — /m\ and v = \ — 5. 

In the result we recognize a special instance of 
the Lauricella functions |l^ defined by 

(33) 



E 



fei 



(ci)fci • ■ • (c„)fc„ fci! • • • fc„! 



fci,...,fe„=0 

where {a)k — r(a + k)/T{a). The defining multi- 
ple series converges for 

. . + < 1 • (34) 



The individual series above converge for mi -f 
™2 + ^/Ip^I < rn-i . Collecting powers in p^ , how- 
ever, the total sum converges due to analyticity 
up to the next singularity on the physical sheet 
given by the threshold condition 

Ip^I < (toi + TO2 + TOa)^ , (35) 

provided that the coefficients themselves do exist, 
which is the case for 



mi + TO2 < "^3 



(36) 



A Lauricella function in the arguments Zi for 
i — 1, . . . ,n can be analytically continued to a 
sum of two Lauricella functions in the arguments 
Xi = Zi/ Zn for i = 1 , . . . , n — 1 and x„ = 1 /z„ by 
the following relation 



F^''\a, 6; ci, . . . , 



(37) 



fiFl^\a,l+a-Cn;ci, . . ., c„^i,l-6+a; xi, ...,„) 
+f2Fl^\b,l + b-Cn;ci, . . ., c„_i,l-a-l-&; xi, ...,„) . 
where 

_ r(c„)r(6-a) _^ 

r(5)r(c„-a)^ ' 
_ r(c„)r(a~6) 

r(a)r(c„ - b) 

Applied to tab. (|l|) this yields, as one coefficient 
vanishes, a total of seven transformed Lauricella 
functions, which is also given in the tab. (^). 
Now, this expression is valid for 

Ip^I > (toi +m2 + m3f . (39) 

One may wonder what the relation is between 
the large p^ expansion of tab. (||) and that given 
in jist . In the latter approach the various terms 
in the expansion are obtained from the expan- 
sion of subgraphs. The subgraphs are obtained 
by distributing the momentum p over the prop- 
agators in all possible ways. In the case of the 
London transport diagram one has the following 
subgraphs: the diagram itself, the three diagrams 
where one internal line is removed and the three 
diagrams where two internal lines have been re- 
moved. 

Following the analysis of |jl^ one can easily find 
the first term of each of the contributing series. 
For the subgraph representing the whole diagram 
the first term in the series should be the mass- 
less diagram. This series then corresponds to the 
last term in tab. @). The series which originates 
from the subgraph where two lines have been re- 
moved, e.g. 1 and 2, starts with the product of 
two massive tadpoles. They contribute a factor 
(rn^ml) which can be identified with the third 
term in tab. (H). The remaining subgraphs are 
obtained by removing one internal line, e.g. line 
3. This yields a series starting with a massive tad- 
pole proportional to (ml)". This is the sixth term 
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Table 1 

The small result 

{zj'z^ T^{-v) F§\l,l + v-l + v,l + u,l + v-zi,Z2,z:i) 
-z'i r2(-zy) 1 - J/; 1 + z/, 1 - I., 1 + zi, 22,23) 

-zl r2(-zy) 1 - 1 - 1 + i^, 1 + t/; zi, 22,23) 

-Y(y)V{^v)V{\ - 2;/) f^^^ (1 - 2t^, 1 - 1 - ^, 1 - 1 + t/; 21, 22, ^3)} 

Table 2 

The large result 

ri23(p^m?,mi,mi) = - [j^y^ {-p')x 
{{-xiYi-xsYT^-iy) F^^\l, l-u-l + u,l-u,l + u- xi,X2. X3) 
+ {-X2Y{~X:iYT^{^v) F^^\l, \^v;\-v,\^v,\^v- x^,X2.x^) 
+ [-xry(^-X2YT'^{-v) F^^\l, + + .ti, X2, X3) 

+ (_.^^)-£!M£l_^lJIi^ - 1 - 2;.; 1 + i^, 1 - 1 - i.; xi, X2, X3) 

rV)ii_gIi^ 4'^ (1 - 1 - 2;.; 1 - J., 1 + J., 1 - Xi, ^2, X3) 
rV)ii_gIi^ 4'^ (1 - 1 - 2;.; 1 - J., 1 - J., 1 + xi, X2, ^3) 
+ 4^(1^3^,l-2t^;l-;.,l-;.,l-^;xi,X2,X3)} ■ 



in tab. (||). Thus the seven series in tab. (|^) can 
be related directly to the seven subgraphs which 
are required for the method of fl^ . 

For the special case p"^ = 0, F^ reduces to the 
Appell function F^ and we recover the result in 
the literature for the vacuum diagram. On 
the other hand taking mi — 7712 — only one 
of the Lauricella functions remains and becomes 
a Gauss hypergeometric function 2-F1 giving the 
result 

ri23(p2; 0, 0, m2) = m2 / __ \ (40) 
r(l-(5)r((5-l)r(2(5- 1)2^1(2^-1, 5; 2-5; z) , 

which is the same as (|2^). 

In the following the general D dimensional ex- 
pression for the 2-loop London transport diagram 
will be expanded in (5 = (4 — Z?)/2 = 1 — for 
small |p^|. 



The following combination |jT^ of the general 
massive case with massless cases is chosen in such 
a way that the infinite parts cancel 

rp / 2 2 2 2\ 

Ti23{p'^;ml,ml,ml) ~ Ti23(p^; m?, 0, TO3) 
-ri23(p'; 0, ml mj) + Ti23(p'; 0, 0, ml) . (41) 

It is this combination which will be calculated in 
two independent ways. 

An analytic form is obtained by expansion of 
the Lauricella functions and their coefficients in i5, 
where the first and the second logarithmic deriva- 
tives of the F-function occur at integer arguments 

n ^ 

i^in + l) - -7 + 51 fc' 
fc=i 

n ^ 

V''(n-f l) = C(2)-Ep- (42) 
fc=i 
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The 1/5^ and 1/5 terms indeed drop out in the 
result and a finite combination of various multiple 
series remains. A good check is provided by the 
cancellation of 7 in this finite expression. 
For small i.e. the region \p^\ < m|, one 
finds 

rr r 2 2 2 2\ / 2 

(i-2)!(t-l)! 



E 



^^_^(m - l)!rn!(n - l)!n!(fc + 1) 



!fc! 



-m n k 



\^{i){t) + i'it-l) - VM - + + log(zi)} 
X {xP{t) + il){t-l) - i^in) - i:{n+l) + log(z2)} 

(43) 

with t = m + n + k. 

4.3. The integral T1234 

With the dispersion method described above 
we derive the small result for ri234 in D di- 
mensions. The discontinuity AT1234 is a sum of 
a two and a three-particle cut, which we denote 

by AT;^^234 AT^^jL' respectively. The two- 
particle cut is given by 

ATi(2) / 2 2 2 2 2\ 

= ABo(p2;TO?,m2)Bo(m2;m2,m^). (44) 

Inserting this result in the dispersion integral 
gives 



1 

27ri 

(mi+m2) 
„2. ™2 



1 



z — 

2\ D 1-^2. „,2 2\ 



= Bo{p ;mi,m2)Bo(™i;"i3'"^4)- 

The discontinuity AT^^L is related to the three 
particle-cut given in ( |27|) but with an additional 
propagator and therefore we get the dispersion 
relation 



n(3) 
-1234 



(p^m|) = -(4V) 



2n2i5 



r2(i - 6) 



(45) 



X yd. 



A 2 °{s,m.-^,mj) 



(5- 



A2 *(s,2,m^) 



ie) J z 

(\/i+m2) 



l-<5 



After performing the z integration we obtain: 



-(A-npL 



2,, r(i-(5)r((5~i) 
r(2-25) 



„<5-l 



rA2 ''(5,7713,7714) 



\47r^2 



F4(l,2-^;2-(5,2-^;^,^) 
s s 



^4(1, 5; 2-5, 5;^,^) 
s s 



(46) 



Expanding (s — 777^) ^ in m1/s and performing 
the integration over s one gets 



l234(P^;"^i) 



r(i -(5)r(i + (5) / 777 



47r/i^ 



-25 



^ Z2"(l - Z3)'^4" 

7n,n,A; J— 

r(i + 777 + r7)r(i + 5 + 777 + ?7 + A: + ;) 
r(2 - 5 + 777)777177!?! 

i_5r(2 - (5 + 777 77)r(2 + m + n + k + l) 
^ T{2-5 + 77)r(2(777, -I- 77 -I- fc) -j- 4 -t- 

r((5 -I- 777 + 77)r(25 + 7r7 + 77 + fc + ?) 



X z. 



r{d + n)r(2(777 + 77 + /c) + 2 + 2(5 + 



(47) 



where Zi = mf/ml with i =^ 1,2,3 and Z4 = 
p^/ml. Note that the contribution from the 
three-particle cut is written in terms of multiple 
series which are not Lauricella functions anymore 
but they belong to a special class of generalized 
hypergeometric functions. To our knowledge they 
have not been studied in the mathematical liter- 
ature. 

With these results ri234 becomes 



71234 (p^;m.^ 



L234(P^"^i) 
2.^2 ^2\u (^2.^2 ^2\ 



(48) 



Since the analytic continuation of the above 
generalized hypergeometric functions is not 
known to us we use the Mellin-Barnes representa- 
tion technique to derive the large p^ expansion. In 
this method a massive propagator is represented 
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in the following way: 



1 1 



ds ■ 



-TO ) 



2\s 



(fc2 - to2)" r{a) 27ri J (fc2)"+« 

— -ioo 

xr(-s)r(a + s), (49) 

where the integration contour in the s plane must 
separate the series of poles of r(— s) on the right 
from the series of poles of r(s + a) on the left. 
In the expression for T1234 we apply (|49| ) with 
a = 1 to all propagators, thereby relating the 
general massive case to the massless one, but 
with the propagators raised to arbitrary powers 
1+Si, I-I-S2, 1+33,1+54. The corresponding result is 
well-known, see e.g.Q. Closing the integration 
contours in a way that the convergence is guaran- 
teed we get seven quartic series each of which can 
again be identified with the subgraph analysis for 
large expansions in |18|. 

The general small result for ri234 can be ex- 
panded in 5. The following combination of the 
general massive case with a massless case is cho- 
sen in such a way that the infinite parts cancel 
|19| 

Ti234:Nip^;ml,ml,ml,ml) = (50) 
Ti234{p'^;ml,ml, TO3, TO4)-Ti234(p^; ml, ml, 0, 0). 

The result can be found in pO| . 

As an example of a special case of these general 
results we choose TO2 = TO3 — and mi — TO4 — 
TO. Eq. ( ^ gives 

r((5)r(2,5)r(l ~5) ( " "^"^ 



1234 ^ (l-S){l-2S) 

X2Fi{2S,S;2~S;x). (51) 

To this the contribution from the two-particle cut 
BQ{p^;m'^ ,0)Bo{m'^;0,m^) should be added 



^ 1234 ^./^ 



r2(i + ^) / 



-25 



<52(l-(5)(l-25) \AT:fi^^ 
X2Fiil,6;2-d;x). (52) 

Expanding in S we get 

T'i234(p^;m-^,0,0, TO^) = 

1 1 
2^ ^ 2^ 



- 2Lrn + 2 ( - — - 1 In (1 - x) 



19 1 

y ~ 2 



C(2) + Ll - 5L„ 



-2L,J ^-^]ln{l-x) + 5(- ^)ln(l-x) 



1 — a;\ 2 / 1 

In [1 — X) Li2{x). 

X I X 



(53) 



We conclude this section with some remarks on 
the master diagram 112345 (p^) in arbitrary num- 
ber of dimensions. This case turns out to pose 
problems. For instance when one wants to apply 
the Mellin-Barnes representation, one needs the 
massless master diagram with arbitrary powers 
for the propagators. This expression is not (yet) 
available in the literature, which is related to the 
fact that the structure of the master diagram is 
not anymore of a self-energy insertion type. 

5. ANALYTIC APPROACHES AND EL- 
LIPTIC INTEGRALS 

In this section we inspect the imaginary parts of 
the London transport diagram T123 and of T'1234. 
It turns out that they can be calculated in four 
dimensions in terms of complete elliptic integrals. 
These are well known functions and thus the re- 
sults are of analytic interest. Furthermore fast 
and precise algorithms for the calculation of the 
elliptic integrals are available. Therefore the re- 
sults provide also an efhcient way to calculate the 
imaginary parts numerically. 

5.1. The imaginary part of the London 
transport diagram 

As can be seen from ( p^ ) the imaginary part of 
T123 is finite in four dimensions and reads with a 
factorization of the Kallen functions 

Im{T,23ip^;m^)) = :^A^l23(p^m2) 

TT f dt r- - 

= J yv(^-a;i)(i-2;2) 

X ^{X3~t)ixi~t), (54) 

Xl = {mi — TO2)2, X2 — (toi + TO2)2, 
X3 = {p- TO3)2, X4, = ip + m^if, 
Xl < X2 < X3 < X4. 



with 
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The integration limits are zeros of the square 
roots, and thus (^J) leads to complete elliptic in- 
tegrals, defined by 
1 

dt 



K{x) = 



E(x) = 



TT 11 2n 



dtjl - xH^) 



(55) 



Il{c,x) = 



TT „ . 1 1 ^^ 



dt 



(56) 



(1 - Ct2)^(l-t2)(l_ 2.2^2) 



|Fi(i;l,l;l;c,x2), 



(57) 



with the Gauss hypcrgeometric function 2F1 and 
the AppeU function Fi ^^^V^. Reduction of 
(^) to the Legendre normal form of the ellip- 
tic integrals |23| , P4| by decomposition into partial 
fractions and partial integration yields after some 
algebra 

Im{Ti23{p'^;ml,ml,ml)) = 

ciK (k) + C2E (k) + can I K 



(58) 



' (toi -1712^ q- + 
+ C4II I , K 



{mi + 7712)2 q 
y.Q {p^ — (mi + 7772 + 7773)^) 



Cl ~ 4777l7n,2 



:|(p + 7773)^ - 7n3P + 777 1 7772 I , 



C2 



C3 



C4 = 



877717772 
877117772 (p^ — "^3)^ 



|(7ni + 7772)(p^ +»773) 

—2m\m\ — 27n3p^|, 



2 94 

K = — 



-9-+ 



with variables q±± corresponding to the physical 
and unphysical thresholds 



q±± [p ± "73)^ - (7711 ± 7712)^ 



(59) 



This result is valid for all values of m\,m\ and 
7M3. In special cases it leads to simpler formulae. 
For equal masses one gets 

/7r7(Ti23(p^; 777^, 777^, 777^)) 
= -^\/{p - m){p + 3777) 



{p — m){p'^ + 3777^) 



E{k) - 47772pK(7t)| 



(p + mY(p — 3m) 



with 



[p — m)^{p + 3m) ' 



(60) 
(61) 



involving only complete elliptic integrals of the 
first and second kind, i.e. 2-^1 Gauss' hypcrge- 
ometric functions. If at least one mass is zero, 
Im{Ti23) reduces to logarithms. 

5.2. The imaginary part of T1234 

The two-particle cut contribution to the dis- 
continuity of ri234 was given in (^J), 

AT^(2) 2 2 2 2 2\ 
^^1234 "72, 7773, 7774) 

^ ABo{p'';mlml)Bo{ml;mlml). (62) 

As a product of a one-loop self-energy integral 
and a one- loop self-energy discontinuity it is com- 
posed of elementary functions and gets a real part 
for 

(7773 -I- 7)74)^ < 777i and (7771 -I- 7M2)^ < P^ ■ (63) 

The three particle cut contribution looks very 
similar to that of the London transport diagram. 
Only one more propagator l/{t — m^) has to be 
added in (p4). The calculation yields 



ATi(3) / 2 2 2 2 2\ 

^-^1234(P im^,m2,m^,m^) 
_ 2ni 



q— 



ciK (k) + C2E (k) + csn I ^— , K 

p^ I \q 

( (7773 - 7774)^ q- + 

-C4II I TTT ,K 



+C5n 



\ (7773 -I- 7774)2 q — 
777i - (7773 - m^ij^ q^+ 
m\ — {ms + 7774)2 q__ 
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xe(p2 



Ci = 4TO3m4 



(64) 



q — 

9++ 



-9 — , 



C3 



8m3m4(p^ 



TOT 



TOS 



ml) 



Ci = 



C5 



8m3TO4(p^ — m|)^ 



mi^q++q — 
8TO3m4A(p^, TO^, m^) 



K — , 

with 

q±± (P ± "^2)^ - ("^3 ± ^4)^. 

In the case (p3) the characteristic c of the last 
IT- function in (p3) is greater than 1, 



(65) 



TO4)^ g_ 



ml — (to3 + m^y 



> 1, 



(66) 



which requires an analytic continuation of that 
function. A comprehensive discussion of the an- 
alytic properties of the elliptic integrals can be 
found in p5| . The ie-prescription in ( |64[ ) en- 

sures that ATI234 gets the correct real part, given 
through 

Im (n(c — ie, k)) 



= — (n(c — ie, k) — n(c + ie, k)) 
2i 



(c- l)(c- k2) 



(67) 



.(2) 



This contribution cancels the real part of Ar£234- 
Consequently Ari234 is always purely imaginary. 

Numerical checks show the agreement of the 
results of (Bq) for /m(Ti23) and of 



/m(Ti234) - Y^ (^^1234 
with previously published tables |16|Jl 



A ^(3) 
1234 



6. ONE-DIMENSIONAL 
REPRESENTATIONS 



(68) 



INTEGRAL 



6.1. A general approach to two-loop inte- 
grals containing a self-energy subloop 

An alternative method to the series expan- 
sion of the two-loop scalar diagrams consists in 



the derivation of one-dimensional integral rep- 
resentations. These are built up from one-loop 
self-energy functions Bq coming from the self- 
energy subloop and the remaining one-loop in- 
tegral. They can be derived by using a dispersion 
representation of the Bq function. 

A two-loop diagram with only three-vertices 




mN+3 



where k is the momentum flowing through the 
self-energy insertion, can in a first step be reduced 
to simpler diagrams by a decomposition into par- 
tial fractions 



1 



1 



1 



X 



fc2 - 



'Af+3 



'^+2 " 
1 



''N+3 



fc2 



'jV+2 



fc2 



''JV+3 



This results for the diagram in 



1 



'■N+2 



''N+3 



X ( Ti,„M+2{Pi,m\, . 
-Ti,,,N+2{p%;rn\, 



2 2 \ 



(69) 



The difference has to be replaced by a derivative 
if TO^_(_2 — nT-^+3- 

Insertion of the dispersion representation for 
the self-energy subloop leads to 



y d-°fcSo(fc^;TO^,m^^J 



[k +piY - m\ 



1 



1 



"N-l 



''N+2 



16 



27ri 



dsABo{s;mjf,m%^^) 



fc2 - 



''N+2 



{k+piY - m{ 



(70) 



After a further decomposition into partial frac- 
tions, 



-1 



1 



k^ 



s k^ 
1 



'Af+2 



m 



N+2 



''N+2 



(71) 



the fc-integrations and one of the s-integrations 
can be performed yielding 



2 2 \ 



Ti...N+2{Pi', "^1, ■ • • , '"■A'+l7 '"-JV+2 



1 

'2TTi 



ds 



ABois;m%,m%^-^) 



'■N+2 



xTW(p,;m^,, 



S). 



(72) 



r(i) denotes a one-loop N-point function in which 
s enters in the remaining one-dimensional integra- 
tion as a mass variable. 

A diagram with two four-vertices leads to a re- 
sult which is similar to the remaining integration 
in ([72I), 
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Pn~i 



rriN+i \mN 



Ti,„N+iipt;m^ 

OC 

1 



27ri 



ds ABo(s;m^,TOjv+i) 



X / d^k 



1 



fc2 — s (fc + pi)^ — m\ 
1 

^ 

(k+pi + ... +PN-if - rri^N-i 

OC 

= ABo(s;TO^,m^^i) 

so 

xT(^)(pj;mi,...,m^_i,s). 
As an example, we apply ( [72| ) to ri234: 

7'i234 (p^ ; , TO2 , mg , to|) 
= Ba{fnl]ml, ■ml)Bo{p'^;ml, mj) 



(73) 



1 
27ri 



ABo(s;mi,TOl) 



(m3+m4)2 

xBo(p^;s,m2). 



(74) 



An application of ( |73| ) to the London transport 
diagram leads to 



123 



1 

27ri 



00 

J ds ABQ{s;m2,m 



(m2+m3)^ 
„2. „ 2\ 



(75) 



a result which would also follow from (|28|). 

The representations (^5|) and ([tJ) with the sub- 
tractions (^) and ( po| ) provide efficient ways to 
calculate T123N and T1234N in a-U parameter re- 
gions. The results agree numerically with those 
published in 



6.2. Integral representation of the master 
diagram 

In this section a one-dimensional integral rep- 
resentation for the master diagram T12345 (fig. |l|) 
in the general mass case is derived. This repre- 
sentation uses only elementary functions and is 
thus well suited for numerical evaluation. 

We use the dispersion relation. The two- 
particle cut contributions are denoted T'^^"^ for 
the cut through the propagators 1 and 2, and 
r(2'') for the cut 4-5. The Cutkosky rules yield 



Ar(2'^)(p2. 



JTii , Too, Too, TO^ 



= ABo(p^\mi,ml 
X (Co(p^ 



2222 2\\* 



(76) 
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This shows that the two-particle cut contribu- 
tions are given by the self-energy discontinuity 
and the complex conjugate of the triangle dia- 
gram Co- For a further evaluation of (|76| ) one can 
introduce the dispersion representation of Co, 



r-i I 2 2 2 2 2 2\ 

ACo(<) 



1 

27ri 



dt- 



Pi 



Coa 



(77) 



with a contribution of the discontinuity belonging 
to the normal threshold 



2m 



:l0g- 



(78) 



VA(i,p^,P^)'"°«-& 
a = t{t + 2mf ^ P2 ^ Ps ^ ™i ^ ^2) 
+ ipl-P2){ml-ml), 

b = \J\{t,pl,pl)^\{t,ml,ml), 

and a contribution Coan belonging to the anoma- 
lous threshold, resulting from the leading Landau 
singularity. This anomalous threshold at 

J. 2/ 2 I 2 I 2 I 2 I 2\ 

ti = m3(— TOg 4- P2 + P3 + '^1 + '^2) 

-{pI -ml){pl-ml) 

-^\{pl,ml, ml)^\{pl,ml, ml) (79) 



occurs if 

mipl -f TO2P3 — m\{'mi + 7712) 
— miTO2(TOi -|- 7712) > 0. 

Following the argument outlined in 
show that its contribution is 



(80) 



one can 



Coa 



dt 



2iTi 



(mi+m2)^ 



Pi ^X{t,pl,pl) 



(81) 



and yields logarithms and square roots. In those 
cases, where is situated on the real axis and 
ti > — (toi -I- 7712)^, this contribution can also 
be taken into account by choosing the appropriate 
sheets for the logarithm in ( |7^). 

After insertion of ( |77| ) into (76) the integrations 
in the dispersion integral can be interchanged and 
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rr(2a) / 2 2 2 2 2 2\ 
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(m4+m5)2 

X (So(p2; 777?, 7772)~Bo(t; 7772,777^)) 



1 

27ri 



(mi+m2)^ 



ABo(s;777j, 777|) 
S — p2 — 76 

,2. _2 2 2\ 



X [Coan{s,ml,ml-ml,ml,ml)) . (82) 

To calculate the contribution of the other two- 
particle cut AT(2''\ the masses have to be in- 
terchanged according to r77i 7774, 7772 <-> 7775. 
Furthermore no complex conjugated form of Cq 
has to be inserted in (|76|), resulting in a complex 
conjugated form of BqIJ,; 777,5) ?a:st term 

of (|8^ ) and no complex conjugated form of Coan 
in the second term of (|2|). 

To evaluate the three-particle cut contribu- 
tions, we keep close to Broadhurst's approach 
1^]. The contribution T^^^) refers to the cut 2- 
3-4, T^^f*) refers to the cut 1-3-5. The Cutkosky 
rules yield 



AT(3-)(p2.^2^^2^^2^^2^^2) 

dt 
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"2^i 



ABoip^;t,ml) 



t- 



76 



(m3+m4)2 
„2 _2. 



xACo(t,p ,7773; 7713,7^4, 77^5). 



(83) 



A partial integration in (p3) leads to a form which 
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simplifies the dispersion integral considerably 

on 

AT(-'''^)(s;m2) 



ds 



I 



s — — le 
ds 



{m2+m3+mi) 
X 



s{s — p2) 



3 — 7712 ) 



X 



xR{s,t,ml,'ml,ml,ml) , 



(84) 



where R is for the case of different masses a ratio- 
nal function of s containing only first order poles. 
An interchange of integrations, 



(v's-7712)^ 

ds j dt 

(7712+7713+7714)^ (7713+7774)^ 



/ / 



dt J ds, 

(7773+7124)2 (m2 + Viy 



(85) 



and decomposition of R{s) into partial fractions 
leads to integrations of the type 



00 

/ 



ds 



y/X{s,t,ml) 



S — Si — It 



= Bo{si;t,m2), 



(86) 



i.e. one-loop self-energy integrals. In the final 
expression all ultraviolet divergences of these Bq- 
functions cancel. The result for the three-particle 
cut contribution of the master diagram is then 



r(3a)(p2.. 
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dt log — 2 - 1 + «e 
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X{t,ml,ml) 
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+Sl(Si - (S3 -I- S4) - S3S4)(sir3 - r^) 
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n (si - Si) 

7=2 

S47-1 + r2 , , 

C4 = — , C3 = C4(S3 S4). 



n (S4 - Si) 

7=1 

Some special mass cases have to be consid- 
ered: double poles occur in R{s) if 7772 = or 
if X{m2,m1,m^) = 0. The decomposition into 
partial fractions leads then to modified results in- 
volving functions 



00 

/ 



ds ^/X{s,t, 
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{s-sf- ief 



dBo{si,t,ml) 
dsi 



(88) 



Another modification of the decomposition into 
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partial fractions occurs if = 0, in which case 
the result is considerably simpler. 

Broadhurst has evaluated the master diagram 
for some cases of physical interest [^. They all 
belong to the special cases mentioned above, if the 
symmetries of the master diagram are exploited. 
The formulae (|8^ ) and provide a method to 
calculate the master diagram in the general mass 
case. The occurring elementary one-dimensional 
integrations are numerically stable and the results 
agree with those of Kreimer's method [|l9|. Com- 
pared with the latter the one-dimensional integral 
representation is much faster, especially if high 
accuracy is required. 

7. CONCLUSIONS 

On the long way to complete two-loop calcu- 
lations in the Standard Model this paper focuses 
on the self-energy diagrams. 

When one restricts calculations to self-energies, 
these should be useful as building blocks with 
suitable theoretical properties. In particular, 
their gauge dependence has to be considered. We 
have discussed this problem for light fermionic 
contributions to the two-loop Z self-energy by ap- 
plying two different approaches, the pinch tech- 
nique (PT) and the background-field method 
(BFM). The PT aims on eliminating the gauge 
parameter dependence of Green functions within 
i?j-gauges, while in the BFM Green functions are 
derived from a gauge invariant effective action. 
Whereas the BFM is valid to all orders in pertur- 
bation theory, the PT had so far been restricted 
to the one-loop level. 

In this paper an extension of the PT to the two- 
loop case has been worked out yielding (within 
i?5-gauges) a gauge parameter independent re- 
sult. It involves a rearrangement of contributions 
between different Green functions and between 
one-particle irreducible and reducible diagrams. 
The BFM has been shown to provide a more gen- 
eral framework. Like at one-loop order, it incor- 
porates the PT results as a special case. The 
BFM vertex functions fulfill simple Ward identi- 
ties which are a direct consequence of gauge in- 
variance. This holds for all values of the quantum 
gauge parameter fg. 



For large parts of the actual calculation of the 
two-loop self-energy diagrams we use an algebraic 
approach allowing for a high degree of automati- 
zation. By applying a method for tensor integral 
decomposition, all two-loop self-energies are re- 
duced to a minimal basis of standard scalar inte- 
grals. 

The evaluation of the scalar diagrams is the re- 
maining problem and is in essence more involved 
than in the case of one-loop diagrams. The rea- 
son is that for the general mass case needed in the 
SM functions beyond the usual (poly) logarithms 
are required. Analytic and numerical approaches 
have been discussed. 

We have derived analytic results in terms 
of generalized hypergeometric functions for the 
scalar self-energy integrals with three and four 
propagators and arbitrary masses. This offers 
the possibility to use well-known mathematical 
techniques like analytic continuation, partial dif- 
ferential equations, or contour integral represen- 
tations. The known formulas for integrals with 
vanishing masses are obtained as special cases and 
consequently are unified in one result in D dimen- 
sions. For the general mass case of the master 
diagram a representation in the form of a gener- 
alized hypergeometric function is not yet known. 

When one is interested only in the imaginary 
parts, an alternative analytic result in four di- 
mensions is obtained in terms of complete elliptic 
integrals. Since their properties are well-known 
they are easily accessible for numerical evalua- 
tions. 

Finally, we have derived a one-dimensional in- 
tegral representation for all two-loop diagrams 
containing a self-energy insertion. For the two- 
loop self-energy diagrams treated in this paper, 
the integrand is composed of elementary func- 
tions only and the representation is valid for all 
values of . Also for the master diagram a one- 
dimensional elementary integral representation is 
derived. The main application of these integrals 
is for numerical evaluations, giving a good alter- 
native to the existing two-dimensional integrals. 

Once adequate techniques for the self-energy 
diagrams are available, one could envisage prac- 
tical applications for physics predictions. For the 
electroweak theory the obvious application is to 
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the gauge boson self-energies which play a role 
in the Mw — Mz mass relation and details of 
the Z line shape. For QED the two-loop vacuum 
polarization is known since a long time |^^, but 
the electron two-loop self-energy was never fully 
calculated. 
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